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We have generalized the dynamical mean-field theory to study the doping dependence of the 
crossover from antiferromagnetic to short-range order modelled by an incommensurate spin density 
wave in the Hubbard model. The local selfenergy which includes spin fluctuations gives quasiparticle 
weights and spectral properties in good agreement with quantum Monte Carlo and exact diagonal- 
ization data in two dimensions. The spectra at finite doping are characterized by a Mott-Hubbard 
'gap' accompanied by a pseudogap induced by the local spin order. 

PACS numbers: 71.27.+a, 75.10.-b, 74.72.-h, 79.60.-i. 
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Although an important progress in the present under- 
standing of strongly correlated fermion systems occurred 
recently due to numerical methods, such as quantum 
Monte Carlo (QMC) and exact diagonalization (ED), an 
analytic treatment that maintains local correlations is 
needed to investigate the thermodynamic limit. An at- 
tractive possibility is the limit of large spatial dimension 
(d = oo), when the diagrams in the perturbative expan- 
sion collapse to a single site and the fermion dynamics 
is described by a local selfenergy This allows to map 
lattice models onto quantum impurity models, which can 
then be solved self-consistently in the dynamical mean- 
field theory (DMFT) 0. At large Coulomb interaction 
U, the spectral properties of Mott insulators were repro- 
duced using the DMFT both at half- (2) , and at arbitrary 
filling n ||, but magnetic long-range order (LRO) was 
analyzed at d = 00 only at half- filling (n = 1) [Q. 

The modification of the magnetic order under dop- 
ing 5 = 1 — n away from n = 1 is a central issue 
in the theory of the oxide superconductors. Undoped 
La2Cu04 is a commensurate antiferromagnetic (AF) in- 
sulator, while doping by Sr into La2_ a; Sr 2 ;Cu04 results in 
short-range AF order within incommensurate magnetic 
structures The physical properties of doped CuC-2 
planes of the high-T c superconductors can be success- 
fully described near n = 1 using a two-dimensional (2D) 
single-band Hubbard (or t — J || ) model, and indeed in- 
commensurate magnetic order was found analytically , 
in Hartree-Fock (HF) S, and in slave-boson approxima- 
tion pPTo|]. Motivated by the successful description of 
the transport properties of the Cu02 planes within the 
DMFT pd],p[, we present in this Letter a generalization 
of the DMFT to the magnetically ordered states, using 



the Berk-Schrieffer spin-fluctuation exchange inter- 
action with an effective potential due to particle-particle 
scattering [[l3| , and construct a selfenergy which leads to 
an overall good agreement with the QMC and ED data, 
ft demonstrates a combination of physics arising from 
the Slater picture and the Mott-Hubbard description of 
strongly correlated electron systems. 

In order to simulate the weakened short-range AF or- 
der in the doped systems, we consider an incommen- 
surate spin spiral (SS) structure. The direction of the 
magnetization m at site i, with respect to the global z- 
axis is specified by two spherical angles, Qi = (4>i,9i). 
Thus, we transform the original fermions {aj^.,at,} to 
the fermions quantized with respect to local quantiza- 



tion axis at each site 



rrA 



a| A , where 



-i(&/2)*, e ^5</2)* y and ~ and b are Pauli 



spin matrices. This transforms the Hubbard Hamiltonian 
with nearest-neighbor hopping, f,-j = —t, to 

H= $3*tf4 ^*(^)ft(%)U' c j(T ,+ Uj2n nnil . (1) 

ij,aa' i 

We take the polar angle as site-independent, 0$ = 0, 
and let the azimuthal angle describe a spin spiral with a 
wave- vector Q, <f>i = Q • Rj. By a straightforward diag- 
onalization of the kinetic energy one finds the energies, 



T Q (k) = - £ k _Q (1 + cos6> a z - sin6> a x ) 



- e k+ Q (1 — cos# a z + sin0 a x ) 



(2) 



where e k = — 2t(cosk x + cos fey). We limit ourselves to 
plane spirals [ful, and choose 8 = n/2. 
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The one-particle Green function is described by a 
(2 x 2) matrix Gq(1c, iu> u ) in spin space, where ui v are 
fermionic Matsubara frequencies. We approximate the 
Green function using a local selfenergy 



(3) 



Gq^k, iw„) = tw„ + /i - f Q (k) - - £ Q F (^„ 



The local lattice Green function, Gq(iuj^) — 
N^ 1 J2 k Gci(k,iuj l ;), becomes diagonal due to the par- 
ity of Tq(k). Therefore, all local quantities includ- 
ing the selfenergy S, which consists of the HF part, 
Sq F = U (no s ) with a — —a, and the spin-fluctuation 
(SF) part, Eq^iccv), are diagonal. Using the cavity 
method for a hypercubic lattice at d — oo, we showed 
that the dynamical Weiss field, Gq a {iu>v) : can be com- 
puted from the Dyson equation of the spin-symmetry 
broken Anderson impurity model, 



Gq («£»>„) 1 



Sq F (^) 



(4) 



The last and the most important step is to find an ap- 
propriate expression for selfenergy. This is known to be 
notoriously difficult and the usual approach within the 
iterative perturbation scheme (IPS) ||] is to develop an 
ansatz which reproduces the exact results in certain lim- 
its. In spite of the correct large U limit ||, the absence 
of a reliable form of E hindered further applications of 
the DMFT for the magnetic states in the intermediate 
U range. By considering the Kadanoff-Baym potential 
containing a class of diagrams up to infinite order p5| , 
we constructed the SF part of selfenergy, 

+ U 2 ^ 1 ^2 X<7<t,q(k<v) SqA^ - iuJ n)> ( 5 ) 

with (3 = 1/ksT. Here we approximated S[G] by S[C/ ]. 
Eq. (||) resembles the selfenergy derived by the cou- 
pling of the moving hole to transverse spin-fluctuations 
p6| , but we note that the coupling to longitudinal spin- 
fluctuations gives a significant contribution. The self- 
energy in a magnetic system is calculated using the 
symmetry-broken magnetic HF state. The transverse, 
Xsa = U[x [ sl] 2 /[^ - Uxfa], and longitudinal, x™ = 
XoS/[l — U 2 x$X^ih susceptibility in Eq. (||) are found 
in random phase approximation (RPA) with renormal- 
ized interaction U. Thus, the renormalizcd interaction 
U is not a fitting parameter ]l7| ], but results from the 

IIP), 



screening by particle-particle diagrams 

U = U/[l + U X pp (0)l 

where x w (0) = &~ x E^ QT (^) 0° Ql (- 
screening effect gives the magnetic structure factor |13|] 
and the selfenergy [jl7| calculated from Eq. (j|) in good 



lOJu 



(6) 



This 



agreement with the QMC results, and depends on the 
underlying magnetic order. It is largest in the paramag- 
netic states and vanishes in the Neel state at n = 1 for 
U — > oo, and is thus very important when the magnetic 
phase diagrams are considered [ p"5| |. 

Eqs. (||), (||), (||), and (||) represent a solution for 
the one-particle Green function within the DMFT in the 
IPS. They have been solved self-consistently and the en- 
ergetically stable spiral configuration was found. This 
procedure is further justified by the sum rule [Q, 

^ " ' '"' " U(n^n ol ), (7) 



-]TE CT ^)G CT ^)e^ 0+ 



which is well fulfilled in the present approach with 
U{n -\nQi) ~ U(not)( n oi) p5[ - We also show below that 
the local selfenergy (^) leads to an overall satisfactory 
agreement with the QMC and ED data. 
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FIG. 1. Total densities of states N(lj) as obtained within 
DMFT (solid lines) for 6 = (AF state), 0.125 [(1,1) spiral], 
and 0.25 [(1,0) spiral] with U/t = 8 and T = 0.05t. The 
dashed lines reproduce N(ui) found by ED of a 4 x 4 lattice 
in Ref. [^2|, and the dotted lines show the Fermi energy. 

The system with local selfenergy is not fully 2D, and 
calculations performed at low temperature T = 0.05t 
(T <C J = At 2 /U) converge to states with magnetic LRO, 
and we discuss here the case of U/t = 8. Away from half- 
filling one finds incommensurate states characterized by 
Q = [tt(1 - 2»fc),7r(l - 2r)y)l V* = Vy for SS(1,1) 
and rj x 7^ 0, r\ y = for SS(1,0), the neighboring spins are 
not strictly antiparallel, the kinetic energy ~ t sin 2?:r] a 
is gained in agreement with polaron physics, and the dy- 
namics becomes closer to free electrons. We reproduce 
a typical pattern found in HF and in slave-boson calcu- 
lations P,|To|]: First, the AF state changes to a mixed 
AF-SS(1,1) state in the range of < 5 < 0.116, which in- 
dicates that a domain- wall phase M might be more stable 
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at low doping. Next, a first order transition to a SS(1,0) 
is found at 5 ~ 0.23, and both spirals coexist in the range 
0.203 < S < 0.248. The moments decrease under doping, 
and finally disappear continuosly at S ~ 0.65. 

Unlike in the IPS based on second order perturbation 
theory the AF order [Q = (ir, ir)] survives at n = 1 for 
U > 4t, as found in QMC at d = oo [^o|. In agreement 
with experiment [ pi) , the spectral functions at n = 1, 
A(k,w) = ImG Q:(J(J ' (k - Q/2,uj + ie), have 

the same generic shape as in t - J model §|19). The 
photoemission spectrum consists of a coherent quasipar- 
ticle (QP) state near the gap with a dispersion ~ 2J, 
and an incoherent part of width ~ 7t at lower energies. 
The DMFT reproduces the nesting instability towards 
the AF order at n = 1, but the AF gap is considerably 
reduced. At U/t = 8 the HF gap of 7.14 t is reduced to 
A ~ 4.93 t, and agrees very well with the ED data for the 
4x4 cluster [p2| (sec Fig. |l|). In contrast, the gap reduc- 
tion is overestimated within approaches [ fl6f not treating 
properly local correlations (pi). 
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FIG. 2. Spectral weight of a coherent QP at k = (ir, 0), 
a(n,0) + 0.2 (triangles), and the averaged weight over Bril- 
louin zone, (a(k))sz (circles), as obtained at n — 1 in DMFT 
(full symbols) and in SCBA (empty symbols) of Ref. 0] . 

The QP weight a(k) increases with increasing J/t, in 
agreement with the self-consistent Born approximation 
(SCBA) (and ED data) [|| in the range J/t < 0.5, both 
at the X = (tt,0) point, and for the averaged weight 
(a(k))ez (Fig. ||). The latter first increases with J/t, as 
in the t— J model, but then flattens out around J/t ~ 0.8, 
and starts to saturate at a value lower than one. This 
behavior follows in first place from the absence of the in- 
tersite SFs, but also shows that the t — J model does not 
apply for J/t > 0.8, where the charge excitations to the 
upper Hubbard band (UHB) become important. 

The quality of our approach is demonstrated by 
the overall shapes of the densities of states N(cu) = 
N^ 1 ^ k A(k, to) which are in a very good agreement with 
the ED data not only at half-filling, but also for doped 
systems (Fig. |l|). At finite doping the system becomes 
metallic, and the selfenergy (^) exhibits a Fermi-liquid 
behavior, with the imaginary part ~ (uj — pi) 2 vanish- 
ing at the Fermi energy \i. This result might change, 
however, in a real system due to critical scattering medi- 
ated by spin-waves. The spectral weight is transferred 



from the UHB into the inverse photoemission (IPES) 
part of the lower Hubbard band (LHB), and N(uj) agrees 
very well with the ED for the SS(1,1) ground state with 
rj x = r\ y ~ S at S = 0.125, both for the positions of 
the UHB and LHB, and their intensities. We reproduce 
thereby the increase of the IPES spectral weight above 
2<5, being 0.331 at U/t = 8 Q. A somewhat smaller 
distance between the LHB and the UHB for the SS(1,0) 
at S = 0.25 might suggest that the system properties are 
closer to a disordered local moment state. 
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FIG. 3. Spectral functions A(k,uj) in the main BZ direc- 
tions: T-X, M-X, and T-M in SS(1,1) state at 5 = 0.125 
and T = 0.05t. The spectra along the V — M direction have 
been averaged over the (1,1) and (1,1) spirals. A shadow band 
in M — X direction is indicated by arrows. 



The k-resolved spectral functions A(k, ui) allow to clar- 
ify the generic features of the doped systems: a large 
Mott-Hubbard gap persists between the LHB and UHB 
(with energies ~ — 3.3t, ~ 6.2t in Fig. |[ respectively), 
while a new smaller pseudogap forms between the pho- 
toemission and IPES parts of the LHB (^5|. The latter 
is induced by the underlying spiral magnetic order which 
removes the degeneracy between X and (7r/2, 7t/2) points 
at S = 0.125, and separates two features which origi- 
nate from a single QP peak at n = 1; the one with a 
higher intensity indicates the position of the QP, while 
the other represents a shadow band. These structures 
and the pseudogap are best visible in the X — M di- 
rection. At the X point one finds the QP peak at en- 
ergy ~ — OAAt and a pseudogap which persists along the 
X — M direction, while the spectral weight moves gradu- 
ally from the lower to the upper maximum |26[ . Thus, the 
QP peak at X point tranforms gradually into a shadow 
band of width ~ 2 J when the M point is approached (see 
Fig. |^) . This interpretation is supported by an excellent 
agreement with the ED data of Ref. Q. The QP in- 
tensities found at 5 = 0.125 are somewhat reduced from 
those in the AF state at n = 1, and agree qualitatively 
with ED for the t — J model at the same doping [fill. 
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Although we cannot attempt a quantitative comparison 
within the t — U model, our results agree qualitatively 
with the QP band with dispersion ~ 2 J and pseudogap 



observed in Bi 2 Sr 2 CaCu208+d' PI 
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FIG. 4. Spectral functions ^(k,^) at S = 0.125 and 
T = 0.33t. The conventions are the same as in Fig. 3. 

The dependence of the QPs on the spiral order explains 
why the ED results obtained for relatively small clusters 
H cannot reproduce the experimental observations that 
the QP peak at X point moves above u> = ji with increas- 
ing doping p7[ . Furthermore, the pseudogap gradually 
disappears with increasing temperature, and the spectral 
weight crosses the Fermi energy along the X — M direc- 
tion, as found at T = 0.33i (Fig. ||). This is consistent 
with the recent QMC calculations |28|], suggesting that 
the pseudogap ~ J originates from short-range magnetic 
correlations at T — > 0. The spectral functions at T = 0.5t 
have only broad maxima which correspond to the LHB 
and UHB, and agree remarkably well with the QMC re- 
sults of Bulut et al. (^9|. At 5 = 0.25 the spectra consist 
of incoherent processes even at low T, accompanied by 
new dispersive features ~ t, which demonstrates that the 
pseudogap and the hole-spin correlations are gradually 
lost in the overdoped regime. 

Summarizing, we have developed a simple but very 
powerful analytic method which captures the essential 
changes of the electronic structure, and gives coherent 
QPs by dressing a moving hole with SFs. While the 
large Mott-Hubbard gap remains almost unchanged with 
increasing doping, the QP states split into two features 
separated by a pseudogap. We believe that the electronic 
properties are generic and not changed by domain-wall 
formation ||. It is expected that the present extension 
of DMFT would significantly improve the understand- 
ing of the strongly correlated transition metal oxides, if 
implemented into realistic band structure calculations. 
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